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Superconducting nanowires fabricated via carbon-nanotube-templating can be used to realize 
and study quasi-one-dimensional superconductors. However, measurement of the linear resistance 
of these nanowires have been inconclusive in determining the low-temperature behavior of phase-slip 
fluctuations, both quantal and thermal. Thus, we are motivated to study the nonlinear current- 
voltage characteristics in current-biased nanowires and the stochastic dynamics of superconductive- 
resistive switching, as a way of probing phase-slip events. In particular, we address the question: Can 
a single phase-slip event occurring somewhere along the wire — during which the order-parameter 
fluctuates to zero — induce switching, via the local heating it causes? We explore this and related 
issues by constructing a stochastic model for the time-evolution of the temperature in a nanowire 
whose ends are maintained at a fixed temperature. We derive the corresponding master equation 
as tool for evaluating and analyzing the mean switching time at a given value of current (smaller 
than the de-pairing critical current). The model indicates that although, in general, several phase- 
slip events are necessary to induce switching via a thermal runaway, there is indeed a regime of 
temperatures and currents in which a single event is sufficient. We carry out a detailed comparison 
of the results of the model with experimental measurements of the distribution of switching currents, 
and provide an explanation for the rather counter-intuitive broadening of the distribution width that 
is observed upon lowering the temperature. Moreover, we identify a regime in which the experiments 
are probing individual phase-slip events, and thus offer a way of unearthing and exploring the physics 
of nanoscale quantum tunneling of the one-dimensional collective quantum field associated with the 
superconducting order parameter. 



I. INTRODUCTION 

The fundamental process governing the collective phys- 
ical properties of quasi-one-dimensional superconducting 
systems is the phase-slip process exhibited by the ex- 
tended, complex-valued superconducting order parame- 
ter field ^'(z), which depends on the position z along 
the system. In the course of a phase-slip process, the 
field ^{z), undergoes a transition from an initial (typi- 
cally metastable) supercurrent-carrying state 5*1 (z) to a 
final one (z). In settings in which the voltage drop 
between the ends of the system is externally controlled, 
these metastable states are topologically distinct from 
one another: the total changes, / dz d^j{z)/dz, in their 
position-dependent phases from one end of the 

system to the other, differ by 2n, and the supercurrents 
carried by these states differ, too 

In the absence of phase-slip processes, the order param- 
eter field of quasi-one-dimensional superconducting sys- 
tems behaves reversibly, i.e., energy stored as kinetic en- 
ergy associated with supercurrent remains undissipated. 
If, however, phase-slip processes do occur, so can dissi- 
pation, part of the coherent kinetic energy of superflow 
being converted into incoherent motion, i.e., heat. For 
example, in a voltage-controlled setting, current-reducing 
phase slips occur with a higher frequency than current- 
increasing ones, leading to energy and current dissipation 
and the notion of an intrinsic resistance, as elucidated 
by Little 0] and Langer and Ambegaokar Q. Similarly, 
in a current-controlled setting [1], the preferred sense of 



phase-slip processes yields an average voltage consistent 
with a positive Joule-heating power. This is the sense 
in which phase-slip processes control the collective prop- 
erties of quasi-one-dimensional superconducting systems: 
they constitute the building blocks via which one can un- 
derstand properties such as dissipation. 

In principle, transitions in the state of the order- 
parameter field can behave predominantly either classi- 
cally or quantally, depending on the temperature of the 
system. In the classical regime, being metastable, the 
states are local minima of the classical free energy, and 
the transitions between such states constitute thermal 
fiuctuations of the order-parameter field over the Arrhe- 
nius energy barriers that separate these states. The study 
of the rates at which such transitions occur, and their 
implications for collective charge-transport through su- 
perconducting nanowires, was initiated by Little 0], and 
developed in detail, shortly thereafter, by Langer and 
Ambegaokar ^] and McCumber and Halperin 5] . In the 
quantal regime, the transitions between the states are 
quantum tunneling events, in which the entire, extended, 
order-parameter field passes from one metastable state 
to another through a classically forbidden region of field 
configurations [aQ- 

In this Paper we develop a theory of the kinetics of 
phase slips of the superconducting order-parameter field 
in settings in which the heat liberated or absorbed dur- 
ing these processes is not instantaneously dissipated but, 
rather, leads to alterations in the local temperature of 
the quasi-one-dimensional system and a resulting fiow of 
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heat, which feed back to influence the phase-sUp kinet- 
ics. As discussed by Tinkham and co-workers [8j], this 
feedback leads to a switching bistability of the system 
involving a pair of mesoscopic states: an essentially su- 
perconducting, low- voltage state and a more highly resis- 
tive, high-voltage state. The rarity of phase slips in the 
essentially superconducting state mean that very little 
Joule heating takes place, which favors the persistence of 
this state. However, the energy liberated by concentrated 
bursts of phase slips can Joule-heat the system enough to 
weaken the superconductivity, which enhances the like- 
lihood of phase slips, and — via this feedback loop — lead 
to the essential destruction of the superconductivity and 
the maintenance of the more strongly Joule- heated, more 
highly resistive state. 

Recent advances in sample preparation techniques 
have made possible the fabrication and exploration of 
extraordinarily narrow nanowires. These wires can be 
so fine, say 10 nm wide or even less, that they bring 
within reach experimental conditions in which heat liber- 
ated during small numbers of phase-slip events — and per- 
haps even a single one — can have a dramatic impact on 
the state of the wire, triggering the switching transition 
from the superconducting to the resistive state. Thus, 
by monitoring the voltage between the ends of the wire 
one should be able to observe macroscopic, or at least 
mesoscopic, consequences of an individual phase slip, 
and hence investigate the properties of these nanoscale 
building-blocks of the collective behavior of quasi-one- 
dimensional superconductivity. For example, one should 
be able to ascertain the rate at which phase slips occur, 
and its dependence on temperature, applied current, wire 
geometry and materials parameters. One should also be 
able to compare such rates with those suggested by theo- 
retical pictures in which the phase slips proceed primar- 
ily via thermal activation over an energy barrier or via 
quantum tunneling through one. Hence, one should be 
able to move beyond the nanoscience that observes the 
structure of nanomaterials or single-particle phenomena 
within them, and progress towards a nanoscience of col- 
lective processes. 

The theory developed in this Paper aims to take a 
step beyond Ref. [1| by considering the stochastic as- 
pects of the phase-slip processes occurring in quasi-one- 
dimensional superconducting systems, i.e., by allowing 
for sequences of phase-slips that occur at random in- 
tervals of time, and exploring the consequences of this 
stochasticity for the states exhibited by the system. Our 
main focus will be on the implications of this underlying 
stochasticity for the rate at which quasi-one-dimensional 
superconducting systems undergo switching transitions 
from the essentially superconducting state to the more 
highly resistive state, as a function of the temperature 
and the current at which the system is maintained. Ex- 
periments are commonly done in a mode in which the 
current is not maintained at a fixed value but is, rather, 
repeatedly ramped up at some fixed rate of increase, the 
current at which the switching transition occurs being 



monitored, so as to produce a distribution of switch- 
ing currents, which depends on the temperature and the 
ramping rate of the current [H, [3 ■ 

Motivation for work reported in this Paper comes from 
the kinds of experimental investigations of superconduc- 
tivity in nanowires touched upon in the previous para- 
graph, and the concomitant need for a road map to guide 
experimental investigations towards regimes of current 
and temperature in which small numbers of phase-slip 
events, or even single such events, induce switching tran- 
sitions of essentially superconducting wires to a highly 
resistive state. Experiments performed in this regime 
should provide access to the temperature- and current- 
dependence of the rate at which individual phase-slip 
events occur. A brief account of this work was reported 
in Ref. UM- 

The Paper is structured as follows: In Sec. [Til we de- 
scribe the switching-current experiments on hysteretic 
superconducting nanowires along with details of our 
physical picture of superconducting-resistive switching. 
We construct a stochastic model of the dynamics of 
the temperature in superconducting nanowires wires in 
Sec. lIIIi and explore its basic properties in Sec.|TVl Next, 
we develop a formalism to address the the statistics of 
switching events in Sec. |Vl and use it to compute the 
switching rate as a function of temperature and bias cur- 
rent in Sec. IVIl and, in turn, compare this rate to exper- 
iments in Sec. I VIII Finally, we present some concluding 
remarks in Sec. IVIIII 



II. PHYSICAL SCENARIO FOR SWITCHING 
IN CURRENT-BIASED NANOWIRES 

The ultranarrow wires that we consider in this Paper 
were fabricated using molecular templating [l^ . By 
using a solution containing long molecules such as car- 
bon nanotubes or DNA one can create a configuration 
in which a nanotube traverses a trench, so as to form 
a bridge-like structure. One can then deposit a layer of 
superconductor, such as MoGe or Nb, on top so that the 
nanotube provides scaffolding on which to form a super- 
conducting nanostructure. In effect, one can thus fabri- 
cate a set-up in which a free-standing superconducting 
nanowire is connected at both of its ends to supercon- 
ducting leads, as shown in Fig. [T^. The diameter of the 
resulting nanowire can be made smaller than the coher- 
ence length of the superconductor, and the length of the 
wire sufficiently greater than the coherence length so that 
the nanowire provides a realization of a quasi-one dimen- 
sional superconductor in which the superconducting fluc- 
tuations are effectively one dimensional. Through careful 
control, the wires produced via molecular templating can 
be made amorphous and quite homogeneous. The result- 
ing superconductor ends up being in the dirty limit (i.e., 
the electron mean-free-path is smaller than both the co- 
herence length and the penetration depth). 

Upon lowering the temperature, the resistance of the 
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nanostructure (i.e., the leads and the nanowire) exhibits 
two drops: a sharp drop as the leads become super- 
conducting, and a second, much smoother transition, 
corresponding to the onset of superconductivity in the 
nanowire itself. The broad resistive transition of the 
nanowire can be understood in terms of the occurrence 
of thermally activated phase-slip fluctuations, and can be 
quantitatively fit in terms of the Langer-Ambegaokar- 
McCumber-Halperin (LAMH) theory 0, [H, Q by using 
the transition temperature and the coherence length as 
fitting parameters. However, the behavior of the resis- 
tance at very low temperatures is not unambiguously 
established, either theoretically or experimentally. On 
the one hand, time-dependent Ginzburg-Landau theory, 
which forms the basis of the LAMH calculation, is not 
strictly applicable in this regime of temperatures and, 
in addition, phase slip processes involving quantum tun- 
neling rather than thermal barrier crossing arc expected 
to become important in this regime. And on the other 
hand, the value of the resistance can fall below the noise 
fioor of the experiment. 

To overcome the difficulty in probing superconductiv- 
ity in nanowires at low temperatures associated with the 
smallness of the linear resistance, we focus on experi- 
ments involving high bias-currents so that they lie be- 
yond the linear-response regime. In these experiments, 
the current through the nanowire is ramped up and down 
in time, via a triangular or sinusoidal modulation pro- 
tocol. As the current is ramped up, the state of the 
wire switches from superconductive to resistive (i.e., nor- 
mal), doing so at a value of the current that is smaller 
than the de-pairing (i.e., equilibrium) critical current; 
and on ramping the current down, the state gets re- 
trapped into a superconductive state, but at a value of 
current smaller than the current at which switching oc- 
curred. Hysteretic behavior such as this, reflecting the 
underlying bistability of the superconducting nanowire 
over a temperature-dependent interval of currents, was 
first reported in Ref. [a]. The experiments addressed 
in the present paper jl2l . ITsj go a step further, in that 
they repeatedly ramp the current up and then down, 
for thousands of cycles at each of a chosen set of tem- 
peratures, and thus generate thousands of values of the 
switching and re-trapping currents at each of these tem- 
peratures. These experiments find that the distribution 
of re-trapping currents is very narrow and does not signif- 
icantly change with temperature. In contrast, the distri- 
bution of switching currents is relatively broad, the mean 
and the width of the distribution changing as the bath 
temperature — which is set by the leads — is varied. The 
fact that even at a fixed temperature and current-sweep 
protocol the switching current is statistically distributed 
and does not have a sharp value is a reflection of (and 
therefore a window on) the collective dynamics of the su- 
perconducting condensate in the nanowire. The conden- 
sate is seen to be a fluctuating entity, evolving stochasti- 
cally in time and, at random instants, undergoing phase 
slip events. The goal of the present work is to understand 



the behavior of these distributions of switching currents, 
and thereby gain insight into the low-temperature rates 
at which thermal and quantum phase-slip fluctuations oc- 
cur. We now proceed to motivate the physical mechanism 
for switching, and thus set the stage for the remainder of 
this Paper. 



Distributions of switching currents were first studied in 
the context of Josephson junctions, in work by by Fulton 
and Dunkleberger [l^. In particular, these researchers 
found that the width of the distribution decreased^ as the 
temperature was reduced. As will be discussed in detail 
in Sec. lVIll in the experiments on nanowires that we are 
considering p^.[l3|. this width is found to increase as the 
temperature is reduced. A second important difference 
is that the Josephson junctions that show hysteresis are 
under-damped systems, whereas nanowires are expected 
to be over-damped, as argued also in Ref. [1], i.e., a single 
phase-slip event by itself is sufficient to cause switching 
in hysteretic Josephson junctions but not in supercon- 
ducting nanowires. Experimentally, the observation of 
voltage tails [13, Hlli [3 ' i-^-, small but non-zero voltages 
across current-biased nanowires in the superconducting 
state, verifies the occurrence of multiple phase slips prior 
to the switching event, and indicates that the wire is 
in the over-damped regime. The main consequence of 
these arguments is that while in the experiments reported 
by Fulton and Dunkleberger the rate of switching is es- 
sentially given by the rate F at which individual phase- 
slips occur, in the case of nanowires the switching rate is 
generically found to be smaller than F. Tinkham and co- 
workers have proposed a physical mechanism that ac- 
counts for the fact that hysteresis in observed, in spite of 
the over-damped dynamics of the wire. According to this 
mechanism, the phase-slip fluctuations are resistive but, 
because of the over-damping, they are not by themselves 
capable of causing switching to a resistive (i.e., normal) 
state. However, the resistance coming from the phase- 
slip fluctuations is associated with Joule heating. If this 
heating is not overcome sufficiently rapidly (e.g., by con- 
ductive cooling) then it has the effect of reducing the 
de-pairing current, ultimately to below the applied cur- 
rent, thus causing switching to the highly resistive state. 
Given that the wire is free-standing, heat can leave the 
wire only by conducting it to the bath provided by the 
leads to which each of its ends is connected. In the fol- 
lowing sections, we discuss in more detail how bistability 
and hysteresis come about within the framework of this 
physical picture. Along the way, it should become clear 
how, within this picture, switching induced by multiple 
phase-slips can essentially be interpreted in terms of the 
number of phase-slips needed to cause a kind of thermal 
runaway instability of the superconducting state. 
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III. BUILDING A MODEL FOR HEATING BY 
PHASE SLIPS 

The goal of this section is to construct a theoretical 
model of the stochastic dynamics that leads to the switch- 
ing of current-biased nanowires from the superconductive 
to the resistive state. We begin by reviewing the theory of 
the steady-state thermal hysteresis as set out in Ref. Q . 
We continue by replacing the steady-state heating of the 
wire by heating via discrete stochastic phase-slips. The 
main result of this section is a Langevin-type stochastic 
differential equation that describes the dynamics of the 
temperature within the wire. 



A. Thermal hysteresis mechanism 

The thermal mechanism for hysteresis in supercon- 
ducting nanowires was originally proposed by Tinkham 
et al. [8|. The qualitative idea of this mechanism, along 
with its relevance to experiments on superconducting 
nanowires, was discussed in the previous section. In the 
present subsection we set up the quantitative description 
of the mechanism by giving a brief account of their work, 
and thus set the stage for our stochastic extension of it. 
Their description rests on the premise that the temper- 
ature of the wire is controlled by a competition between 
(i) Joule heating, and (ii) cooling via the conduction of 
heat to the baths (leads). If Q{x) is the temperature at 
position X along the wire of length L and cross-sectional 
area A, then the power-per unit length dissipated due to 
Joule heating at a bias current / is taken to be 
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l'R{Q{x),l) 



(1) 



where the function i?(8') is to be understood as the 
resistance of an entire wire held at a uniform temper- 
ature 8(x) = Q' . On the other hand, as the wire is 
suspended in vacuum, the heat is almost exclusively dis- 
sipated through its conduction from the wire to the su- 
perconducting leads that are held at a temperature Tj, 
and which play the role of thermal bath. The heating 
and cooling of the wire is described by the corresponding 
static heat conduction equation: 



Qs 



,{x) = -d^[K,{Q)dMx) 



(2) 

Ks{e)dl&, (3) 



where Ks{Q) is the thermal conductivity of the wire (The 
first term on the right-hand side of Eq. ([3]) was absent in 
Rcf . [8] ) . This equation is supplemented by the boundary 
conditions 6(±L/2) = Tb at the wire ends, x — ±L/2, 
and we solve it numerically via the corresponding dis- 
cretized difference equation. 

It was found in Ref. [S] that Eq. ^ yields two solutions 
for a certain range of / and Tb. The nonlinear dependence 



of the resistance R on temperature, which is characteris- 
tic of a superconducting nanowire, is at the root of this 
bistability. This bistability in turn furnishes the mech- 
anism for the thermal hysteresis in the I-V characteris- 
tic; the two solutions correspond to the superconducting 
(cold solution) and the resistive (hot solution) branches 
of the hysteresis loop. To obtain the hysteresis loop at a 
given bath temperature Tb, one begins by solving Eq. ^ 
to obtain Q{x) at a bias current sufficiently low such that 
the equation yields only one solution. Next, by using 
the solution Q{x) from the previous step to initialize the 
equation solver for the next bias current step, the locally 
stable solution of Eq. ([2]) is traced out as a function of / 
by tuning the bias current first up and then down. The 
I-V loop is thus traced out by calculating the voltage. 



V 



L/2 



L/2 



dxIR{e{x),l), 



(4) 



at each step. 

The numerical analysis of Eq. ([2]) requires a knowledge 
of i?(6,/) and Ks[&), which serve as input functions 
for the theory. A discussion of these input functions and 
other parameters is given in Appendix!^ In Ref. Q, the 
linear-response resistance measured at Tb = was used 
for i?(0, /). However, i?(8, /) depends also on the value 
of the bias current /. Moreover, we find that by incor- 
porating deviations in R{Q,I) from the linear-response 
regime, we are also able to obtain a better fit with the 
experiments considered in this Paper (see Ref. [T3[ and 
Appendix |^ . 



B. Heating by discrete phase slip events: 
Derivation of Langevin equation 

In the previous subsection, we described the static the- 
ory of thermal hysteresis as was discussed in Ref. @ in 
the context of experiments on MoGe nanowires. Let us 
now go one step further and include dynamics by consid- 
ering the time-dependent heat diffusion equation: 

Cv(e) 9*9(2;, <) = a, [K,ie)d^eix,t)] + Qsourcc, (5) 

where the specific heat Cv(0) enters as an additional 
input function. This differential equation can be derived 
in the standard way, by using the continuity equation. 



^ ■ JQ + dtQ ^ Qsour. 

for the heat current, 

jQ = -if,(e)ve, 

together with the energy density. 



Q 



Cy{e')de'. 



(6) 



(7) 



(8) 



However, as long as we assume that the wire is heated 
by the source term given by Eq. ([1]), the dynamic for- 
mulation turns out to be inadequate for our purposes, as 
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a. Experimental setup 



bath 



bath 



b. Phase slip 
Re \|/ 
Im 

c. Simplified model 



d. Temperature profile 
T 



FIG. 1: Schematic showing a. the configuration of the free- 
standing nanowire supported by the superconducting leads, 
which act as a thermal bath; b. the attenuation of the order 
parameter tp{x) in the core of a phase slip; c. the simplified 
model of the wire: phase-slips occur exclusively in the central 
segment while the end segments carry the heat produced by 
phase-slips to the thermal baths; d. the temperature profile 
of the wire in the simplified model, with uniform temperature 
in the central segment and spatially varying temperature in 
the end segments. 



should become clear from our analysis and results. Such 
a source term assumes that the wire is being continually 
heated locally as a result of its resistivity R{'d{x), I)A/ L 
at any given position x along the wire. Is this assump- 
tion of continual Joule heating correct? To answer this 
question we need to deconstruct the resistance and get 
to its root. In Sec. |TT] we have dwelt upon phase-slip 
fluctuations in detail. There, we have emphasized the es- 
sential point that it is the resistive phase-slip fluctuations 
that are responsible for the characteristic resistance of a 
quasi-one dimensional wire. Thus, one should consider 
the Joule heating as being caused by individual, discrete 
phase-slip events. 

Let us then explicitly consider discrete phase-slip 
events (labeled by i) that take place one at a time at 
random instants of time ti , and are centered at random 
spatial locations Xi. By using the Josephson relation 



2eV I-kV 



dt 



(9) 



relating the voltage pulse V{t) to the rate of change of the 
end-to-end phase difference across the wire (j), we arrive 
at the work done on the wire by a phase slip, viz.. 



where $o = h/2e is the superconducting flux quantum. 
Hence, a single phase slip (or anti-phase slip), which cor- 
responds to a decrease (or increase) of cj) by 2tt, will heat 
(or cool) the wire by a "quantum" of energy Wps- By us- 
ing this result we can now write down a time-dependent 
stochastic source term: 



x^)5{t-ti), (11) 



where F{x — Xi) is a spatial form factor, of unit weight, 
representing the relative spatial distribution of heat pro- 
duced by the i^^ phase-slip event, and ai — ±1 for phase 
(anti-phase) slips. The probability per unit time T± for 
anti-phase (phase) slips to take place depends on the lo- 
cal temperature Q{x, t) and the current /. 

Now, instead of using the continual Joule-heating 
source term, Eq. Ilj, let us use the source term given 
by Eq. pT|) . Instead of being a deterministic differen- 
tial equation, the heat diffusion equation ([5]) becomes a 
stochastic differential equation for T{x,t). We thus have 
a Langevin equation with stochasticity, in one space and 
one time dimension, with a "noise" term that is charac- 
teristic of a jump process. 

Let us pause to understand the connection between the 
two source terms. By using the Josephson relation ([9|. 
we can express the resistance as 



i?(e,/) 
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<i>or(e,/) 



(12) 



and use it to rewrite the continual Joule-heating source 
term, Eq. ([T]), as 



AL ' 



(13) 



where F = F_ — F+ is the net phase-slip rate for the entire 
wire. Let us assume that a phase slip only affects its local 
neighborhood, i.e., F{x — Xi) ~ 5{x — Xi). Then, if we 
take the continuous time limit of Eq. (jlip by assuming 
that the phase-slips are very frequent and that Q-ps — > 0, 
the two source terms would indeed become equivalent (as 
can also be seen formally by taking the limit 'I'o ^ 0). 

We now make a brief remark about the switching cur- 
rent. The static theory of hysteresis that was discussed 
in the previous subsection has a single, well-defined value 
of the switching current, which corresponds to the value 
of the bias current at which the low-temperature (su- 
perconductive) solution becomes unstable. On the other 
hand, we see from the theory discussed in the present 
subsection that the randomness in Xi and ti generates a 
stochasticity in the switching process. The full implica- 
tions of the stochastic dynamical theory will be discussed 
in the following sections. 



C. Simplified model: Reduced Langevin equation 



= / dtIV = / 



4 = *o/, 



(10) 



In principle, one can proceed to study the physics of 
the stochastic switching dynamics of a current-biased 
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nanowire by using the dynamics of Eq. ([5]) together with 
stochastic source Eq. (ITT|) . both derived in the previous 
subsection. In practice, however, it is not easy to solve 
the full Langevin equation with both spatial and tempo- 
ral randomness. In this subsection we derive a simpli- 
fied model, and argue that it is capable of capturing the 
physics essential for our purposes. 

We concentrate on wires that are in the dirty limit, 
for which the mean free path is much shorter than the 
coherence length, which is shorter than the charge im- 
balance length required for carrier thermalization, which 
itself is somewhat shorter than the nanowire length L. In 
addition to restrictions on length-scales, we assume that 
the time for a phase-slip (~ tgl) and the quasi-particle 
thermalization time te are both smaller than the wire 
cooling time, i.e., the time it takes the heat deposited in 
the middle of the wire by a phase-slip to diffuse out of 
the wire. 

We will make a series of simplifications as follows: 

1. Due to the presence of the superconducting leads 
at two ends, as well as edge effects, it is more likely 
that the phase-slip fluctuations in the wire are cen- 
tered away from the wire edges. We thus assume 
that the source term is restricted to a region near 
the center of the wire. 

2. We assume that the heating takes place within a 
central segment of length I, to which a uniform tem- 
perature T is assigned. Note that the total length L 
may be allowed to differ slightly from the geomet- 
ric length of the wire, in order to compensate for 
the temperature gradients in the lead at the wire 
attachment points. 

3. We assume that heat is conducted away through 
the end segments, each of which are of length (L — 
Z)/2. As an additional simplification, we ignore the 
heat capacity of these end segments. 

4. To simplify the problem further, we make use of 
the fact that the probability per unit time (i.e., 
the characteristic rate) r+ for an anti-phase slip 
to take place is much smaller than the rate r_ for 
a phase slip to take place, and we thus ignore the 
process of cooling by anti-phase slips. To account 
indirectly for their presence, we use a reduced rate 
r = r_ — r+ instead of r_ . This ensures that the 
discrete expression for Q correctly reduces to the 
continual Joule-heating expression. 

With the simplified model defined above, the descrip- 
tion of superconducting nanowires reduces to a stochas- 
tic ordinary differential equation for the time-evolution 
of the temperature T of the central segment: 

^ = -a(T, Tb) (T - Tb) + rKT, /) ^ 6{t - U). (14) 

i 

This equation can be regarded as a spatially reduced 
Langevin-type equation which is a counterpart to the 



full, spatially dependent Langevin-type equation Q and, 
correspondingly, has a spatially reduced version of the 
full source term (jlip . The second term on the RHS of 
Eq. corresponds to (stochastic) heating by phase 

slips, and the first to (deterministic) cooling as a result of 
the conduction of heat from the central segment to the ex- 
ternal bath via the two end-segments. The temperature- 
dependent cooling rate a(T, Tb) is obtained by compar- 
ing the heat currents through the end segments to the 
thermal mass of the central segment, where the heat 
currents through the end segments are found by solv- 
ing the heat equation in the end segments, subject to 
the boundary conditions that T(0) — T{L) — Tb and 
T(-^) = T{L - ^) = T. The cooling rate may be 
expressed in terms of the integral 

"(^'^-)" /(L-Oa(T)T-Tb X/^'^-(^^^- ^''^ 

If Ti and Tf are temperatures of the central segment be- 
fore and after a phase slip then we can express the tem- 
perature 'impulse' due to a phase slip, i.e., Tf — Tj = 
77(Ti, /) = rj{Tf, I), as function of either Ti or Tf (depend- 
ing on the context) by using 

A I c,{r)dr = Qp,. (16) 

JTi 

To summarize, in this section we have derived a sim- 
plified model that is described by the reduced Langevin 
equation (fT4| . The central assumption that we used to 
build this simplified model is that phase-slips predomi- 
nantly occur in the center of the nanowire, or at least 
that their exact spatial locations are unimportant. This 
assumption is appropriate for shorter nanowires, in which 
we do not have several distinct locations along the wire 
at which a switching event may nucleate. Specifically, 
if the wire length does not greatly exceed the charge im- 
balance length (which itself is assumed to be much larger 
than the coherence length) then, independent of where a 
phase-slip occurs, the temperature profile in the wire af- 
ter the phase slip will be similar, and we are able to apply 
our simplified time-only model. 



IV. BASIC PROPERTIES OF THE SIMPLIFIED 
MODEL: BISTABILITY AND SWITCHING 

The goal of this section is to explain the basic prop- 
erties of the stochastic model that was formulated in 
the previous section and is described by the reduced 
Langevin equation, Eq. (|14p . First, we show explicitly 
that the competition between the heating and cooling 
terms leads to the emergence of bistability. Next, we 
show that the stochastic character of the heating term 
can lead to switching between the two metastable states. 
Then, to characterize this switching we relate the lifetime 
of the superconducting state to the Mean First Passage 
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FIG. 2: Effective potentials for various bias currents (0.175, 
0.195, 0.215, 0.235, 0.255 ^A) at fixed lead temperature Tb = 
1.2 K [13] • The inset shows the details of the local maximum 
of the effective potential, and corresponds to an enlargement 
of the region indicated by the box in the main plot. 




Temperature 

FIG. 3: Schematic depiction of the effective potential for the 
bistable case, showing the labeling scheme for the various tem- 
peratures: the bath temperature Tb, effective-potential min- 
ima corresponding to the superconducting state Tac and the 
resistive state Trs, the local maximum of the effective poten- 
tial Tap, and the Mean First Passage Time point T*. 



Time (MFPT) for the temperature T in the central seg- 
ment to exceed a certain critical value T*. 

We begin by reminding the reader that the theory of 
Ref. @, described by Eq. ^ with the source term given 
by Eq. ([T]), is the static continual-heating version of the 
theory described by Eq. ([5]) with the source term given by 
Eq. (fTTj) . The connection is made evident via Eq. ([13]). 
As has been shown at the end of Subsection IIII Al for 
certain values of Tb and / the continual-heating theory 
describes a bistable system having a low-temperature su- 
perconducting branch and a high-temperature resistive 
branch. The fluctuating theory described by Eqs. ^ and 
(llip . in turn, allows for processes that move the system 
between the two metastable states. 

Analogously, there is a simplified continual-heating 
theory associated with the spatially-reduced theory de- 
scribed by the reduced Langevin equation (|14p . To obtain 



the continual- heating version of Eq. ((14)) , in analogy to 
Eq. (|13p . we replace the term carrying the sum over the 
delta functions by the phase slip rate r(r, /). That is, 
if we ignore the discreteness of phase-slips, the dynamics 
of the temperature is described by 



dt 



= -a{T, Tb) (T - Tb) + 77(T, /) r(T, /). (17) 



We may think of this equation as describing the motion 
of an over-damped "particle" of position T{t) at time 
moving in the fictitious potential U{T), i.e., 



dT 
'dt 



dUjT) 
dT ■ 



(18) 



In Fig. [2] we have plotted the form of this potential for 
several different values of the current at Tb = 1.2 K for 
parameters corresponding to a typical nanowire. At low 
values of the bias current, the fictitious potential U{T) 
has only one minimum, which corresponds to the super- 
conducting, low-temperature state with T T^. As the 
bias current is increased, a second minimum, correspond- 
ing to the resistive state, develops at higher tempera- 
tures, and the system becomes bistable. For the bistable 
regime we label the temperatures of the two minima of 
the effective potential by Tgc and Tj-g, and the tempera- 
ture of the local maximum that separates them by Tgp, as 
depicted schematically in Fig. [3l Further increase of the 
bias current results in the high-temperature minimum 
gradually becoming deeper and the low-temperature min- 
imum shallower. Eventually, the low-temperature min- 
imum disappears and only the high-temperature mini- 
mum remains. 

Consider a system biased such that it is bistable in 
the continual heating description. If the fluctuations due 
to the discreteness of phase-slips are weak, i.e., the tem- 
perature rise caused by an individual phase slip is small 
compared to the temperature difference between the two 
metastable minima, then the system highly likely to re- 
main in whichever of the two minima it started in. How- 
ever, very rarely, the intrinsic fluctuations in the times 
between phase-slips will drive the system from one local 
minimum to the other. 

A picture of a switching event can be constructed by 
analyzing the real-time dynamics of the Langevin equa- 
tion. A typical trace of the temperature of the central 
segment as a function of time T(t), evolving according to 
Eq. (fn)) . is depicted in Fig. H) To obtain this trace we 
integrate Eq. (|14p forward in time, numerically, starting 
with the initial condition T{t = 0) = Tb. Phase slips 
correspond to sharp rises of T(t), and cooling to grad- 
ual declines of T{t). The two minima, Tgc and Tjs, as 
well as the local maximum Tgp of the fictitious potential 
U{T), are indicated by the red lines. From the trace, it 
can be seen that the temperature of the system starts by 
spending a long time in the vicinity of Tgc, until a burst 
of phase slips pushes it "over" Tgp and the temperature 
quickly progresses towards the vicinity of Tig. 
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FIG. 4: Trace of a typical central segment temperature tra- 
jectory T{t), with a switching event at f ~ 20 ps [TtI . Over- 
laid on the trace are red lines showing the temperatures T^c 
(lowest), Tap (middle), and Trs (highest) corresponding to the 
superconducting minimum, the saddle-point, and the resis- 
tive minimum of the effective potential. The inset shows a 
blow-up of the boxed region of the main figure. 



The fundamental quantity of interest is the mean 
switching time, i.e., the average time required for the 
wire to switch from being superconductive to resistive. 
Assuming that the entire wire has temperature T = Tb 
when the current I is turned on at time t — 0, we define 
a switching event as the first time at which T, the tem- 
perature of the central segment of the wire, exceeds the 
temperature T*, where Tgp < T* < T^-s. With this defini- 
tion, the mean switching time corresponds to the Mean 
First Passage Time T(Tb T*) to go from the bath 
temperature Tb to the temperature T*. In the case of 
weak fluctuations, the problem is indeed a barrier cross- 
ing problem, i.e., the system spends a long time in the 
vicinity of the starting temperature Tb ~ Tgc, until a 
burst of phase slips propels it over the barrier at Tgc. Af- 
ter Tsc is exceeded, the system moves quickly (compared 
to the barrier-crossing time) towards Tig. Therefore, the 
mean switching time will only have a weak dependence on 
T*, provided T* is significantly higher than Tgp, as this 
is the temperature range in which the system is moving 
relatively quickly towards the high-temperature steady 
state. 



V. FORMALISM FOR ADDRESSING 
SWITCHING DYNAMICS: MEAN FIRST 
PASSAGE TIME 

In the previous section, we have shown that the Mean 
First Passage Time (MFPT) for the temperature T in the 
central segment to exceed a critical value T* can be used 
to characterize the switching from the superconducting 
to the resistive state. In this section, we develop the 
tools for computing the MFPT in two steps. In the first 
step, we derive the Master Equation associated with the 



Langevin equation [141 In the second step, we use the 
Master Equation to obtain a delay-differential equation 
directly for the MFPT r(T) [= t(T T*)] as a function 
of the initial temperature T. 

A. Master Equation 

We now derive the Master Equation associated with 
the Langevin equation ^4]) . The Master Equation is a 
delay-differential equation that describes the evolution 
of the probability density P{T, t) for the central segment 
of the nanowire to have temperature T at time t. Note 
that by the term "delay" we actually mean a delay in 
temperature. That is, the evolution of P(T,t) is non- 
local in temperature but local in time. Before delving 
into the derivation, we quote the result: 

dtP{T, t) ^dr [{T - Tb) a{T) P{T, t)] - T{T) P{T, t) 
+ T{T-^{T))P{T-7,{T),t)x 

X {l-dT?j{T)), (19) 

where the first (i.e., the transport) term corresponds to 
the effect of cooling, and the last two terms correspond 
to the effects of heating. We remind the reader that, 
rj{T) is defined as the temperature change due to a phase 
slip, such that after the phase slip the central segment 
temperature becomes T. 

The probability distribution P(T, t + dt) at tempera- 
ture T and time t + dt is related to P{T, t) at an earlier 
time t by the two effects: (a) cooling, and (b) heating, 
as summarized in Fig. [S] To understand these effects, 
we begin by discretizing the temperature interval into 
equally wide slices Tj of width dT, indexed by i. 

To understand the effect of cooling, we consider the 
change in the probability of finding the system in the i-th 
temperature slice, P{Ti,t + dt) — P{Ti,t). Due to cooling, 
some of the probability in slice i + 1 will move into the 
i-th slice at a rate (T^+i — Ti,) a{Ti^i) / dT . Concurrently, 
some of the probability slice i will move into slice i — 1 
at a rate {Ti — Tb)a{Ti)/dT. These two processes are 
indicated by the blue arrows in Fig. [S] By adding these 
two rates, wc find 

^tP{T^,t) = ^[{T,+,-n)a{T,+,)P{T,+,,t) 

-{T,~n)a{T,)P{T„t)]. (20) 

To find the continuum version of this equation we take 
the limit dT 0, thus arriving at 

dtP{T,t) ^ dT[{T -n)a{T)P{T,t)], 

where we have identified the definition of the T derivative 
when taking this limit on the right hand side of Eq. (|20p . 

On the other hand, heating is caused by discrete 
phase-slips, which occur at a temperature- and current- 
dependent rate r(T, J). Heating decreases the probabil- 
ity in the i-th slice, P(Ti, t), at the rate r(Ti, /) P{Ti, t), 
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FIG. 5: Schematic diagram showing the probabihty P{Tj,t) to find the system in the temperature interval {Tj, Tj + AT) 
indicated by gray blocks. Processes involved in the rate of change of the probability in the i — th interval (green block) are 
represented by blue arrows (cooling) and red arrows (heating by phase-slips). The boundary conditions at the ends of the 
temperature interval (Tb, T*) for computing the Mean First Passage Time are indicated in purple. 



as the probability is boosted to higher temperatures by 
phase-slips. On the other hand, P(Tj,t) increases, as 
probabihty from lower temperatures, T « Ti — rj(Ti), gets 
boosted to Ti due to heating. These two heating pro- 
cesses are indicated by the red arrows in Fig. O To com- 
pare these two rates, we must take care of the fact that 
the boost is temperature dependent, and thus there may 
be some 'stretching' of the corresponding temperature 
intervals before and after the boost. Consider the tem- 
perature interval {Ti, Ti + dT) after the boost. What is 
the corresponding temperature interval before the boost? 
Using the 'un-boost function' rj, we find the interval to 
be 

(T, - ^(T,), T, + dT- ?y(T, + dT)) . (21) 

Therefore, the width of the interval before the boost is 
approximately (1 — v{T^i))d'T and not dT . Thus, after 
taking the continuum limit we find that the probability 
density P{T) increases at the rate 

(1 - dr viT)) f{T - ?iiT)) P{T - ?,iT), t). 

Combining the rates of in- and out- flux of the probability- 
density due to cooling and heating effects, we obtain the 
master equation given as Eq. 

B. Mean First Passage Time 

The fundamental quantity that we want to compute 
is the Mean First Passage Time (MFPT) T{Th T*). 
In this subsection, we use the master equation ([T9|) . as 
a starting point for obtaining a delay-differential equa- 
tion directly for r(T) = t{T T ). We proceed usmg a 
straightforward generalization of the standard procedure 
to systems with jump processes (i.e.. those that have mas- 
ter equations with delay terms) [H, [l^ . 

We begin by supplementing the Master Equation with 
the boundary conditions appropriate for computing the 
MFPT. When computing the MFPT, we want to re- 
move any element from the ensemble once its tempera- 
ture reaches one of the boundaries of the interval. Thus, 
we would typically impose absorbing wall boundary con- 
ditions on both sides of the interval at Tb and T*, i.e., 
P(Tb,t) = P{T*,t = 0). However, because we have a 



jump process, which transfers probability-density from 
lower temperatures to higher temperatures, we must in- 
stead impose the absorbing segment boundary condition 

P{T,t) = 0, forT>r*, (22) 

beyond the upper end of the interval, so as to capture 
systems in which the temperature gets boosted beyond 
the upper end of the interval by the jump process. The 
dynamics described by Eq. pop is also peculiar in another 
way. If we choose Tb as the temperature at the lower 
end of the interval then there are no processes that can 
cause a passage through the lower end of the interval, as 
can be seen from the Langevin equation p^ . Thus, at 
the lower end of the interval, instead of the absorbing 
wall boundary condition, we must impose the no flux 
boundary condition 

aT,P(Tb,i) =0. (23) 

These boundary conditions, ((22|) and ((23|) . are indicated 
(in purple) in Fig.O We note that the probability density 
may be discontinuous at the upper boundary, due to the 
jump process. 

With the boundary conditions in hand, we consider the 
integrated probability density function 

G{T',t)= f dTP{T,t\T\0). (24) 

Here, we have generalized from the probability P{T, t) 
to the conditional probability P{T, t\T', 0) of finding the 
system to be at temperature T at time t, subject to the 
boundary conditions, if it started out at temperature T' 
at time 0. Then the function G(T', t) measures the prob- 
ability that a system that started out at temperature 
T' has never left the temperature interval (Tb, T*) while 
the time t has passed. In particular, the rate of first pas- 
sages out of the interval (Tb, T*) at time t is given by 
-dtG{T',t). Therefore, the MFPT t(T') for a system 
that starts at temperature T' is given by 

t(T')= / dtt{-dtG{T',t)) (25) 
Jo 

= / dtG{T',t), (26) 
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where the surface term resulting from an integration by 
parts is assumed to be zero, as ah the "particles" are 
assumed to be able to leave the interval in the long-time 
limit. 



Next, we obtain a differential equation for t{x) by 
appropriately integrating the backwards-in-time master 
equation, i.e., the equation 



dt'P{T, t\T', t') = (T' - Tb) a(r') dT'P{T, t\T', t') + r(T') [P(r, t\T' , t') - P{T, t\T' + r]{T'), t')] 



(27) 



Taking advantage of the fact that the present stochastic 
process is homogeneous in time, we transfer the time- 
derivative on the left hand side of Eq. to the t vari- 
able from the t' variable: 



dt'P{T,t\T',t') 



-dt'P{T,t 
-dtP{T,t~ 



-t'\T',0) 
t'\T',0). 



(28) 
(29) 



By substituting Eq. p9|) into Eq. (|27|) and integrating 
both sides with respect to T' over the interval (Tb, T*), 
we arrive at an equation for G{T',t): 



dtG{T', t) = (Tb - T') a{T') dT'G{T', t) 

- r(T') [G(T', t) - G(r' + r/(T') 



t)] 



(30) 
(31) 



Finally, we integrate over all times to obtain an explicit 
delay-differential equation for the MFPT: 

{n-r)a{r)dT'T{T') (32) 

-r(T') [r(r')-r(T' + 7y(T'))]=-l, 

where we have used Eq. (^5)) to identify t{T') on the left 
hand side, and the assumption that P{T, t\T\ 0) tends to 
zero in the long-time limit on the right hand side. 

The delay differential equation (|32p . together with the 
boundary conditions (l22|) and ([SS]) . may be conveniently 
solved numerically by using the shooting method. The 
key to this method lies in taking advantage of the fact 
that in the nonlocal term r(T') r(r' -t- ?7(r')), the factor 
r]{T') is always positive. Therefore, by integrating from 
high temperatures to low temperatures, we can always 
look up the value of the non local term from the region 
where integration has already been carried out. We im- 
plement the shooting procedure as follows: (1) pick a 
value for t(T*); (2) shoot towards lower temperatures to 
obtain dT^T(T\,); and (3) adjust t{T*) until the boundary 
condition dr^, t{T\^) = is satisfied. 



VI. SWITCHING BEHAVIOR OF 
SUPERCONDUCTING NANOWIRES: RESULTS 
AND INTERPRETATIONS 

A. Properties of solutions of the Mean First 
Passage Time equation 

The Mean First Passage Time t(T' T*) for a sys- 
tem that starts out at temperature T' to exceed the tem- 



I=0.2yUA 



I=0.465yuA 




1.2 



2.81.2 



2.8 



I=0.485yuA I=0.545/iA 




I=0.625yuA I=0.725yL/A 




u 



2.81.2 



T(K) 



T(K) 



2.8 



FIG. 6: Numerical solutions of the Mean First Passage Time 
Eq. (I32|l . as functions of the central segment temperature T, 
for bath temperature Tb = 1.2 K and several values of the 
bias current between 0.2 and 0.725 (solid pink) [T^]. For 
comparison, superposed on the MFPT plots are the effective 
potential that appears in Eq. (|18p (dashed red lines). The 
solid blue lines correspond to the zero of the effective poten- 
tial. The hash marks on the temperature axes correspond to 
the sequence of temperatures given by Eq. 



perature T* is described by the delay-differential equa- 
tion ([311), together with the boundary conditions ([^ 
[23)) . The mean switching time corresponds to the MFPT 
T(Tb T*). To understand the solutions of the 

MFPT equations, we remind the reader that T* must 
be chosen to be a temperature sufficiently far above the 
saddle-point temperature Tgp of the effective potential 
(cf. Fig. [3]) that the MFPT has only a weak dependence 
on T* . Solutions of the MFPT equations for several val- 
ues of the bias current are plotted in Fig. [51 along with the 
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corresponding effective potentials. The solutions have 
the following structure: 

• In the temperature interval Tb < ^ < Tgp the 
MFPT is largely independent of the temperature 
T. 

• For T in the vicinity of T^p the MFPT drops 
sharply. 



In the temperature interval Tgp < T < (T* , T^s 
MFPT is again largely independent of T. 



the 



The origin of this structure can be seen in the real- 
time dynamics depicted in Fig. |4l Systems that start at 
temperatures below the barrier temperature Tgp in the 
effective potential must diffuse over it, which is a very 
slow process. The MFPT in the interval Tb < T < Tsp is 
correspondingly large. Furthermore, the MFPT is largely 
independent of the initial temperature, as this interval is 
essentially ergodic, i.e., a system that starts in this in- 
terval typically spends a lot of time exploring the entire 
interval before leaving it. On the other hand, a system 
that start at some temperature above the barrier in the 
effective potential rolls down the potential gradient rela- 
tively quickly before reaching T*. Therefore, in terms of 
temperatures increasing from Tb, the MFPT starts out 
essentially constant over the ergodic interval, and then 
drops sharply as T crosses the barrier in the effective 
potential, before finally flattening out at temperatures 
higher than T^p. 

We note that it is numerically advantageous to set T* 
to be as low as possible in the high-temperature (i.e., 
flat MFPT) regime, so as to avoid instabilities in the 
numerical integration of Eq. 



B. Number of phase slips in a thermal runaway 
train 

In this subsection, we consider the question of exactly 
how many phase-slips it takes to form a runaway train. 
Let us start our discussion by first turning off the (de- 
terministic) cooling term in the stochastic equation (|14p . 
If we now start with an initial temperature Tq, then the 
sequence 



To, To + v{To), To + v{To) + v{To + v{To)), 



(33) 



defines the discrete sequence of values that phase-slips 
would cause T to jump to, as marked on the horizontal 
axes in Fig. [S] for Tq = Tb. The probability per unit time 
r(r) to make a jump changes at each step, and so does 
the size ri{T) of the jump, owing to their explicit depen- 
dence on temperature. On the other hand, if we turn 
off the heating term then we would have a deterministic 
problem in which T would decay at a rate o:{T), from 
its initial value Tq > Tb to the bath temperature Tb. It 
is the competition between the discrete heating and con- 
tinuous cooling that makes for a rather rich stochastic 
problem. 




/(mA) 

FIG. 7: Map of N{Tb, I), i.e., the number of consecutive phase 
slips required to overcome the barrier and induce the wire to 
switch, as a function of bias current / and bath temperature 
Tb- The three solid lines are contours along which (the 
mean switching rate) is 1, 10'^, and lO'' s~^. Similarly, the 
three dashed lines are contours along which the phase slip- 
rate is 1, 10'^, and 10® s~^. For comparison, the de-pairing 
critical current is indicated by the dashed-dotted line. 



The number of tick marks [see sequence p3p] between 
Tb and T (see Fig. [6]) is nothing but the of number 
N(T\y,I) of phase-slip events required to raise the tem- 
perature of the central segment from Tb to T in the 
absence of cooling. Accordingly, iV(Tb,/) also provides 
an estimate of the number of phase-slip events needed 
to overcome the potential barrier, if the time-span of 
these events were insufficient to allow significant cool- 
ing to occur. "Thermal runaway" — heating by rare 
sequences of closely-spaced phase slips that overcome 
the potential barrier — constitutes the mechanism of the 
superconductive-to-resistive switching within our model. 
As the number iV(Tb,/) of phase-slips needed increases, 
the total number of phase-slip events taking place before 
a switching event occurs, and correspondingly the value 
of the mean switching time Tg (Tb, /), may indeed become 
quite large. 

The map of the number iV(Tb, /) of phase slips needed 
over the TTb plane is presented in Fig. [71 The progres- 
sion from smaller to larger iV(Tb,/) corresponds to the 
progression from lighter to darker shading. For clarity 
of presentation we have grouped together a few values of 
Af(Tb, /). We see from the figure that, as the bias current 
is increased, one traverses through regions of decreasing 
A^(Tb, /). It is important to observe that the typical value 
of iV(Tb,/) is only ten or fewer. The smallness of this 
number highlights why going beyond continuous Joule 
heating to the discrete phase-slip model is crucial for our 
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FIG. 8: Logarithm of the mean switching rate, InFsw, 
as a function of current for bath temperatures Tt = 
0.1, 0.3, 0.5...,1.9K. 



analysis. 

Remarkably, there is a region in the /-Tb plane within 
which the occurrence of just one phase-slip is sufficient 
to cause the nanowire to switch from the superconduc- 
tive to the resistive state. We denote this this region 
as the "single phase-slip switching regime" . A switching 
measurement in this range can in fact provide a way of 
detecting and probing an individual phase-slip fluctua- 
tion. 



C. Mean switching time 

Let us begin by considering the single phase-slip 
switching regime identified in the previous subsection. 
In this regime, the value of the mean switching time is 
dictated purely by the probability for a phase-slip event 
to occur at a given bath temperature Tb and by the bias 
current /. The mean switching rate is thus identical to 
the phase-slip rate (which is an input quantity in our 
theory) : 



0.0 0.25 0.5 0.75 1.0 1.25 

/(mA) 

FIG. 9: Switching rate Psw as a function of current for bath 
temperatures Tt = 0.1, 0.3, 0.5...,1.9K. The current sweep 
rate was set to 58/j,A/s. 



contour lines decreases monotonically upon lowering Tb, 
the spacing between the lines can be seen to behave 
non-monotonically. In the following subsection we exam- 
ine a consequence of this non-monotonicity. 



D. Switching-current distribution 

The mean switching time Tg in bistable current-biased 
systems can be measured directly by performing waiting- 
time experiments. Alternatively, 7= can be extracted 
from switching-current statistics |16j . As described in 
Section ini the switching-current distribution can be gen- 
erated via the repeated tracing of the I-V characteristic, 
ramping the current up and down at some sweep rate 



dl 
'di 



±r. 



The sweep-rate-dependent probability P(I,Th;r)dI for 
the event of switching (from the superconductive to the 
resistive branch) to take while the current is in the range 
I to I + dl is explicitly related to the mean switching 
time r= via the relation 



"single phase-slip switching regime." (34) P(Tb I-r)dI— fr^^(Tb /) 



As we move beyond the "single phase-slip switching 
regime," several phase-slip events become necessary for 
switching. The mean switching rate thus begins 
to deviate from F to a value below F, as can be seen 
from Fig. [51 For concreteness, we have assumed that the 
phase-slip fluctuations are thermally activated, and have 
used the form of F as given by the LAMH theory for the 
current-biased case. (The precise expression is provided 
in Appendix \^) The deviation of from F can be 
interpreted in terms of the evolution in the number of 
phase slips A^(Tb,/), which was discussed in the previ- 
ous subsection. To make the link between the two more 
evident, we have color-coded the plots according to the 
value of A^(Tb,/). 

An alternative representation of our results can be ob- 
tained by studying the equal-value contours. A graphical 
representation of the contour lines for a few values of 
and F, chosen in an experimentally accessible range, is 
provided in Fig. [T] Whilst the spacing between the F 



dl 



X 1^1-^ P(Tb,/V) d/'j . (35) 

The term in the first pair of parentheses corresponds to 
the probability for switching to happen within the ramp 
time, whilst the term in the second pair of parentheses 
corresponds to the probability that the wire has not al- 
ready switched before reaching the bias current J. By 
using Eq. ()35|) we obtain the distribution of switching 
currents in superconducting nanowires in terms of the 
theory presented in the present paper. The plots illus- 
trated in Fig. [9] are simply the translation of the mean 
switching-rate plots shown in Fig. [BJ to the switching- 
current distribution for a chosen value of r. 

Upon raising Tb, one would naively expect the dis- 
tribution of the switching currents to become broader 
for a model involving thermally activated phase slips. 
Such broadening of the distribution is indeed obtained 
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up to a crossover temperature scale T"{r) (i.e., the tem- 
perature below which, loosely speaking, switching is in- 
duced by single phase slips). However, on continuing to 
raise Tb, but now through temperatures above r"(r), 
the distribution-width shows a seemingly anomalous de- 
crease. Based on the discussion of the previous subsec- 
tion, this can be understood as a manifestation of the 
now-decreasing spacing between the Tg contour lines in 
the /-Tb space. 

This striking behavior above T"{r) may be understood 
via the following reasoning: the larger the typical number 
of phase-slips in sequences that induce superconductive- 
to-resistive thermal runaway, the smaller the stochas- 
ticity in the switching process and, hence, the sharper 
the distribution of switching currents. This non- 
monotonicity in the temperature dependence of the width 
of the switching-current distribution, along with the ex- 
istence of a regime in which a single phase-slip event can 
be probed, are the two key predictions of our theory. In 
the following section, we proceed to carry out a detailed 
comparison of between this theory and the recently per- 
formed experiments discussed in Sec. |TT1 



4: T^=0.3K 



VII. COMPARISON WITH EXPERIMENTS 



In this section, we compare results from the exper- 
iments described in Section |TT] with predictions of our 
theory presented in Sections IIIHIVll We show that our 
theory is both qualitatively and quantitatively consistent 
with experimental observations. The main implications 
of this comparison are that: (1) the switching-current 
distribution-width does indeed increase as the temper- 
ature is decreased; (2) there is a single phase-slip-to- 
switch regime at low temperatures; and (3) thermally 
activated phase-slips, alone, are insufficient to fit the de- 
pendence of the mean switching time on the bias current 
at low temperatures. This suggests that one should in- 
clude the effects of Quantum Phase-Slips (QPS) 0, Q; 
upon including quantum phase-slips phenomenologically, 
we obtain good fits to the experimental data in the low- 
temperature regime as well. For the purposes of this 
comparison we use the data from a representative su- 
perconducting nanowire; data for more samples may be 
found in Refs. [H, [l^l- This section is structured as fol- 
lows. To establish the validity of the thermal hysteresis 
model, we begin by analyzing the I-V hysteresis loops. 
Next, we qualitatively analyze the experimentally mea- 
sured switching-current distribution. We continue with 
a quantitative analysis of the experimental data on the 
mean switching rate. Finally, we look at the implica- 
tions of the quantitative analysis, including identifying 
the single phase-slip-to-switch regime and the scenarios 
for quantum phase-slips in the low temperature behavior. 
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FIG. 10: Comparison of I-V hysteresis loops computed us- 
ing central segment approximation and the full heat equa- 
tions. For the case of the heat equation the wire was split 
into 15 segments. The steps and upturns in the high volt- 
age branch of the segmented solution correspond to various 
segments becoming superconducting, and would disappear in 
the continuum limit. Heat conductivity, phase-slip rate, and 
geometrical parameters used are listed in Table U 



A. I-V hysteresis loops 

We begin our analysis by comparing the qualitative 
features of the experimentally measured and theoreti- 
cally computed current-voltage characteristics. Experi- 
mentally, it is found that at high temperatures there is 
no hysteresis. As Tb is lowered, a hysteresis loop grad- 
ually opens up. Next, as the temperature is lowered 
even further, the switching current (i.e., the bias current 
at the superconducting-resistive transition) grows gradu- 
ally, whilst the re-trapping current (i.e., the bias current 
at the resistive-superconducting transition) remains al- 
most unchanged. This behavior is consistent with the 
experimental observations and theory of Ref. 8] , where 
it is also qualitatively explained as follows. Switching is 



controlled by the properties of the low-temperature (i.e., 
superconducting-like) solution of Eq. , thus switching 
depends strongly on the temperature of the bath T^. On 
the other hand, re-trapping is largely a property of the 
hotter (i.e., resistive) state, and thus has only a weak 
dependence on Tb. 

Typical I-V curves obtained from the central segment 
model [e.g., stationary solutions of Eq. (fTT])]. as well as 
those obtained from solving the heat equation, are shown 
in Fig. (TU] for several bath temperatures Tb. Following 
Ref. , the solutions of the heat equation were obtained 
from a spatially discretized version of Eq. In both 
cases, the heat conductivity and the phase-slip rate were 
obtained from Eqs. (jAlsp and (|A1[) . respectively. The 
theoretical curves both qualitatively and quantitati vely 
reproduce the features seen in experiments [1, [13, Il3| . 
Wc take a moment to point out that in fitting the ex- 
perimental data it is important to take into account the 
nonlinear dependence of the phase-slip rate on the bias 
current. Finally, we point out that making the central 
segment approximation has little effect on the switching 
current found in the hysteresis loops (for typical wires 
used in experiments). This fact supports the validity of 
the central segment approximation for modeling switch- 
ing phenomena. 



B. Switching current distributions 

In the experiments, every time an I-V characteris- 
tic is measured by sweeping the bias current up and 
down, switching occurs at a distinct value of the bias 
current. By repeatedly measuring the hysteresis loop 
at a fixed bath temperature Tb and current sweep rate 
dl/dt, one can obtain the distribution of switching cur- 
rents P{h^,T^,dI/dt). Typical P{I^^,T^,dI / dt) distri- 
butions, obtained experimentally, are shown in Fig. [TT] 
For completeness, we also show the corresponding theo- 
retical fits, which we shall describe in detail in the next 
subsection. For a given Tb, switching events tend, in gen- 
eral, to occur at lower bias currents than the switching 
current found in the thermodynamic stability analysis of 
Ref. Q. The reason for this premature switching at bias 
currents that are lower than the stability analysis indi- 
cates is, of course, thermally activated barrier crossing in 
the form of a phase-slip bursts. In Fig. [12] we plot the 
mean and the standard deviation of the switching-current 
distributions measured experimentally, as well as those 
obtained from theoretical fits of the simplified model. By 
using the tuning parameters obtained from the fits, we 
also plot the theoretical de-pairing critical current and 
the critical current from the stability analysis of the sim- 
plified model Eq. pT]) . 

As described in the introduction, one would typically 
expect the standard deviation of the switching-current 
distribution to decrease with decreasing temperature, as 
thermal fluctuations become suppressed. Such narrowing 
of the switching-current distribution is expected to con- 
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FIG. 11: Experimental data (red bar-charts) and theoreti- 
cal fits (black lines) for the switching- current distributions 
for various bath temperatures. To make the comparisons be- 
tween the shapes of distributions easier, we also show theo- 
retical curves shifted so that their means coincide with the 
experimental curves (blue dashed lines). The fitting parame- 
ters used are listed in Table |T] 
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FIG. 12: Hierarchy of critical and switching currents. All the- 
oretical curves were produced using the parameters obtained 
from fitting experimental data. The highest current scale 
is the de-pairing critical current Ic{T). Next, is Iu{T), the 
scale for linear instability due to overheating, as described by 
Eq. (|17|) . Finally, comes the actual switching current Jaw(r). 
The error bars on the switching current correspond to the 
10 X the standard deviation of the switching current distribu- 
tion (the scale of the standard deviation was exaggerated to 
make it easier to see; the sweep rate was set to 58 /^A/s). The 
fitting parameters used are listed in Table U 
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FIG. 13: Standard deviation of the switching-current distri- 
bution as a function of temperature. Comparison of experi- 
mental data for a typical sample (circles) with various QPS 
scenarios labeled by T, . The "best fit" curve was obtained by 
using the parameters given in Table HI where the fit was opti- 
mized to simultaneously capture the temperature-dependence 
of mean switching current and the standard deviation of the 
switching current distribution. 



tinue with cooling, until the temperature becomes suffi- 
ciently low such that quantum phase slips are the main 
drivers of the switching, at which point the narrowing 
is expected to come to a halt. Qualitatively, this would 
indeed be the case if switching was always triggered by a 
single phase-slip. However, our theory, predicts that the 
situation is more complicated, because the mean switch- 
ing time, and hence the width of the switching-current 
distribution, is controlled by a competition between the 
phase-slip rate and the number of consecutive phase-slips 
needed to induce switching, as described in the previous 
section. Thus, qualitatively, we expect the opposite be- 
havior at higher temperatures. That is, in the regime 
of thermally activated phase slips and at temperatures 
above the single phase-slip-to-switch regime, the width 
of the switching-current distribution should increase with 
decreasing temperature. This counter-intuitive broaden- 
ing of the switching current distributions with decreasing 
temperature is indeed observed experimentally, as shown 
in Fig. M 



C. Mean switching rate 

As the switching-current distribution depends on the 
bias-current sweep-rate, in order to quantitatively com- 
pare our theory with experimental data, we focus on the 
mean switching rate, which is related to the switching 
current distribution via Eq. (|35p . The experimentally ob- 
tained mean switching rate t^^(/, Tb) for a typical sam- 
ple, along with theoretical fits, are plotted as a function 
of the bias current / for different values of the bath tem- 
perature Tb in Fig.[T4l To help relate the mean switching 
rate to the switching current distribution width, we note 
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FIG. 14: Comparison of theoretical fits (lines) to the exper- 
imental (blue dots) Mean Switching Rates. The top panel 
shows fits with thermally activated phase slips only, whilst the 
bottom panel includes both thermally activated and quantum 
phase slips. The solid red lines correspond to multiple-phase- 
slips-to-switch regime, whilst the dashed blue lines correspond 
to multiple-phase-slips-to-switch regime. The fitting parame- 
ters used are listed in Table |T] 



that for a fixed Tb, the shallower the slope of t^^{I,T\)) 
the wider the corresponds distribution. The two main 
features of the experimental data plotted in Fig. [M] are 
as Tb decreases (1) the mean switching rate decreases 
((/sw ) increases) and (2) the slope of t~^{I) decreases 
(/sw distribution width becomes wider). 

Two different fits to the same set of experimental 
data are shown in Fig. [T31 The fit shown in the top 
panel includes TAPS only, whilst the one in the bot- 
tom panel uses the fitting parameters from the top panel 
but also includes QPS. The fits were obtained using the 
fast switching-rate calculation routine described in Ap- 
pendix [C] The tuning parameters that were obtained 
from the fit are listed in Table [T] and fall into two cate- 
gories. The first category is composed of the geometric 
model parameters, such as the wire length, whilst the 
second category controls the "input functions," i.e., the 
heat capacity, the heat conductivity, and the phase-slip 
rate. The expressions for these input functions are given 
in Appendix \X\ We note that in obtaining these fits we 
verified that the fitting parameters that we used were 
consistent with the high-temperature R{T) data (l3| . 
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The TAPS-only fit (top panel of Fig. [H]) works well 
at temperatures above IK. In this regime the theory is 
able to quantitatively explain the observed rise in mean 
switching current ((/sw)) with decreasing temperature, as 
well as the peculiar increase the Igw distribution width 
with decreasing temperatures. 



D. Single-slip-to-switch regime 

In general, as the temperature is lowered and the bias 
current increased, the wire tends to enter the single-slip- 
to-switch regime, as indicated in Fig. [71 This regime 
roughly corresponds to the region of the (/, T) plane 
where a single phase slip heats up the wire to Tc{I), and 
thus the boundaries of this regime are primarily deter- 
mined by the heat capacity of the wire. Within this 
regime, switching-current statistics correspond directly 
to the phase-slip statistics. 

Theoretical fitting indicates that at temperatures be- 
low ~ 1 K the wire enters the single-phase-slip-to-switch 
regime. This regime is indicated by the switch of the 
theory curves from solid red lines to dashed blue lines in 
Fig. 1141 In the absence of quantum phase slips, in this 
regime the Igw distribution width should follow a more 
conventional behavior, and decreases with temperature. 
This corresponds to the increase in the slope of the mean 
switching rate curves with decreasing temperature in the 
single-slip-to-switch regime (see the top panel of Fig. [T4l) . 

However, experimentally the distribution width seems 
to increase monotonically as the temperature is lowered, 
even in the single-slip-to-switch regime. This behavior 
suggests that there is an excess of phase slips at low tem- 
peratures. 



E. Quantum phase-slip scenarios 

We expect that at low temperatures quantum phase- 
slips will contribute strongly to the switching rate. We 
model the presence of QPS by adding their rate to the 
rate of TAPS, so as to obtain the total phase-slip rate 
which goes into our model: 

rtotai(/, T) = rTAPs(/, T) + rQPs(/, t). 

To model the QPS rate, rQPs(/, T), we replace ksT 
kBin+TiT) in rTAPs(/,r) (see AppendixEI])- Here, 
and Ti are both treated as fitting parameters. Letting 
Ti be nonzero does somewhat improve the quality of our 
fits. 

We can envision several scenarios for the effect of QPS 
on the switching current distributions, and these are 
summarized in Fig. 1131 In the absence of QPS, upon 
lowering the temperatures, once the single-slip-to-switch 
regime is reached the distribution width will start de- 
creasing with temperature. This type of behavior is 
demonstrated by the T, = line in Fig. [T31 However, 



in the presence of QPS, the distribution width is ex- 
pected to saturate at low temperature, with the satura- 
tion value controlled by T». If, upon cooling, the single- 
slip-to-switch regime is reached before the temperature 
reaches T*, we expect the distribution width to first in- 
crease and then decrease before saturating with decreas- 
ing temperature (cf. the = 0.5 K curve in Fig. [T3|) . 
On the other hand, if is reached before the single-slip- 
to-switch regime is reached, we expect the distribution 
width to increase monotonically with decreasing temper- 
ature (cf. the ^ 1.0 K, 1.5 K curves in Fig.[l3l). 

To include QPS in our fitting, we started with param- 
eters obtained by fitting the mean switching rate curves 
at high temperatures (T > IK), as described in the pre- 
vious subsection (i.e. see top panel of Fig. [H)) . Next, we 
optimized T* and Ti to obtain the best possible fit to the 
mean switching rate curves at low temperatures as well. 
The optimal values thus obtained were T, = 0.726 K and 
Ti = 0.4, which corresponds to the fit shown in the bot- 
tom panel of Fig. [14] and the curve labeled "best fit" in 

Fig.m 

To fit the experimental data, we must be able to simul- 
taneously match both the mean and the standard devi- 
ation of the switching current distribution as a function 
of temperature. However, we have not been able to get 
quantitative agreement with both of these quantities, si- 
multaneously, in the low temperature regime. The pa- 
rameter values of ~ 0.726 K and Ti = 0.4 result in 
a good fit of the mean but not the standard deviation 
(see Figs. [131 and [131), whilst the values of T, = 1.5 K and 
Ti — result in a good fit of the standard deviation (see 
Fig. [T51) but not the mean (not shown). 

We conclude this section by noting that, for the 
nanowire that we fitted, our fitting seems to favor the 
QPS scenario where T* is higher than the temperature 
corresponding to the onset of the single-slip-to-switch 
regime. 



VIII. CONCLUDING REMARKS 

In conclusion, we have developed a quantitative theory 
of stochastic switching from the superconducting to the 
normal state in hysteretic superconducting nanowires. 
Our theory describes the dynamics of the heating of 
nanowires by random-in-time phase-slip events, and the 
cooling of the wire by heat conduction into the leads. In 
general, a train of phase slips, sufficiently closely spaced 
in time can cause the nanowire to overheat and switch 
from the low-temperature (superconducting) branch to 
the high-temperature (normal) branch. 

The main achievement of our theory is that it quantita- 
tively describes the unexpected increase in the switching 
current distribution with decreasing temperature, as is 
observed in experiments. 

Our theory also predicts that there is a single-slip-to- 
switch regime at low temperatures. In this regime, a 
single phase slip always triggers a switching event; thus. 
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by studying the switching statistics one has direct access 
to the phase shps statistics. Typically, phase-slip proper- 
ties have been studied by linear-response measurements, 
which are only feasible when the nanowires have a mea- 
surable resistance, i.e., at high temperatures. The single- 
slip-to-switch regime is interesting because it occurs at 
low temperatures, and thus it is a complementary tool 
with which to study the properties of phase slips. 

Finally, the monotonic increase of the switching cur- 
rent distribution width with decreasing temperature, 
even in the singie-slip-to-switch regime, seems to indi- 
cate a severe excess of phase slips over the predictions 
of the Langer-Ambegaokar McCumber-Halperin model 
of thermally activated phase slips. It seems at the very 
least plausible that the quantum tunneling of the super- 
conducting order parameter (i.e., QPS) is the mechanism 
that serves to meet this excess jia. 
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APPENDIX A: INPUT FUNCTIONS AND PARAMETERS 



In this appendix we catalog the models for the phase-slip rate, heat capacity, and thermal conductivity that go into 
the stochastic heat equation 



1. Phase-slip rate 



We begin by considering the phase-slip rate. For the TAPS rate, Ptaps, we have used the LAMH model, including 
the nonlinear current response: 

rTAPs(/,r) -r_,TAPs(/,T) -r+,TAPs(/,r), (Al) 

r±,TAPs(/,r) - f}±(/,T)exp (- ^^^^y^^ ) , (A2) 

where -I- or — indicate whether the phase slip results in current rise or drop, respectively. The phase-slip barriers 
AF± at bias current / and temperature T are given by 



AF^{T)=CiiT) (l^Vl - 3A:2 - 8fc(l - fc^) arctan j > (A3) 



AF+{T) = Ci{T) -\/2VT~~3F + 8fc(l-fc2) 



TT — arctan ^ 

V2k 



(A4) 



C^i(T) = ^^/.(r), (A5) 

where the phase gradient k at current / is the real solution of the equation 

/ 



Ic{T) 



= k{l-k^), (A6) 



and the temperature-dependent critical current Ic{T) |20l | is expressed in terms of the wire length L, the critical 
temperature Tc, the zero-temperature coherence length and the normal-state resistance of the wire i?„ [21|, via 



/ / ^ 2\ 3/2 
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We approximate the prefactor ^±{1, T) in Eq. (jA2p via 



n+{i,T) ^n{T). (Aio) 

In the presence of a bias current /, the "+" phase shps are exponentiaUy more rare than the "-" phase shps. Therefore, 
we keep the current-dependent corrections to the prefactor for the "-" phase shps, but not for the "+" phase shps. 
Thus, we obtain an approximation that works in both the hnear-response regime, where the current correction is 
irrelevant, and in the high-bias regime, where phase shps are rare. We estimate the temperature-dependent 
coherence length and the Ginzburg-Landau relaxation time via 



Thus, we can express the phase-slip rate via the physical parameters L, i?„, ^c- To obtain the quantum phase- 
slip rate, we replace AF/ksT in in Eqs. and Eq. by AF/ksTes where T^s = (T, + TiT) is the effective 
temperature. and Ti are treated as fitting parameters, T* being the low-temperature limiting value of Tcff. 



2. Heat capacity and thermal conductivity 

Unfortunately, we know of no direct experimental data on the heat capacity and thermal conductivity of current- 
carrying superconducting nanowires. The diameter of the wires used in experiments is comparable to ^o- Thus, the 
thermodynamic properties of these wires should lie somewhere between those of a bulk superconductor and a normal 
metal. Therefore, for the purposes of computing the thermodynamic functions, we model the wire as being composed 
of a BCS superconducting wire of cross-sectional area Ai in parallel to a normal-metal wire of cross-sectional area 
A2. The BCS and Fermi liquid expressions for heat capacity [22l | are 

a,FL(e) = l^'^NoklQ, (AM) 



where (3 = l/fcs6, Ft — ^J^+S^Q), fk is the Fermi function, and A(6) is obtained from the BCS gap equation. 
Thus, the total heat capacity of the wire Cy is given by 

- -mTa, • ^^^^^ 

Similarly, the dirty-limit BCS [23l [2^ and Fermi-liquid expressions for thermal conductivity are 

7 A 2/cBe ksO 

^.flO)^^, (A17) 

where D is the diffusion constant (for MoGe D ~ Icm^/s [2^), and Lq = ir^kg/Se^. The total thermal conductivity 
Kg is, correspondingly, given by 

A1+A2 ^ ' 

The fitting parameters describing the heat capacity and thermal conductivity are the cross-sectional areas Ai and 
A2, and Tc of the nanowire. 

APPENDIX B: FITTING PROCEDURE switching rates at high bias currents and low tempera- 



The main goal of the fitting procedure is to fit the 
switching-rate data. However, in addition to the mean 
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type 


parameter name 


symbol 


value 


geometric 


wire length 


L 


110 nm 


geometric 


central segment length 


Li 


110 nm 


geometric 


end segment length 


L2 


27.5 nm 


input function 


transition temperature 




3.872 K 


input function 


zero temperature transition length 


Co 


5 nm 


input function 


QPS effective temperature 


Teff- = T* +TiT 


0.726 K + 0.4T 


input function 


effective superconducting cross-sectional area 




320.4 nm^ 


input function 


effective normal cross-sectional area 


A2 


19.0 nm^ 


input function 


normal state resistance 


Rn 


2666 STi 



TABLE I: Parameters used in the switching rate model. Parameters fall into two categories: geometric parameters and input 
function parameters. The former define the simplified model of the wire, while the latter define the phase slip rate, heat capacity 
and heat conductivity. Not all of these are used as tuning parameters, as J?„ and L can be measured directly, and the fitting 
is only weakly effected by Li and L2. Fitting parameters for a representative experimental sample (corresponding to Figs. 1111 
1121 1131 114p are displayed in the right hand column. 



tures, we also have data on the linear-response resistivity 
in the high-temperature regime. (The linear response re- 
sistivity becomes too small to measure below T ^ 1.9 K 
for our wires.) The fitting is performed in two steps. In 
the first step, we fit the high-temperature linear-response 
data. In the second step, we use the parameter val- 
ues from the first step as a starting point in fitting the 
switching-rate data. Table |T] lists the parameters that go 
into our model; the procedure for determining them is 
explained below. 

We fit the high-temperature linear-response by con- 
ductivity following the usual procedure f21'| . In this pro- 
cedure, L and i?„ are obtained from microscopy and elec- 
trical measurements of the wire resistance above Tc- We 
fit the R{T) data using 

ii; = lim -5^ = lim J $0 rTAPs(/, T) 
1^0 I i-^o I 

to obtain Tc and ^o- 

Next, we use the values of L, _R„, a-nd Tc, obtained 
in the first step, as a starting point in fitting of the mean 
switching-rate data. In this step, we tune Ai, A2, Tc, 
£_o, T* , and Ti simultaneously to obtain the best pos- 
sible fit over the entire current and temperature range. 
During this procedure, we set Li = L and L2 = L/A. 
We find that variation of Li and L2 does not signifi- 
cantly effect the fit, and thus we exclude them from the 
already-extensive list of fitting parameters. 

Finally, we verify that the fitting parameters obtained 
from fitting the mean switching-time data are consistent 
with the high-temperature linear-response data. 



APPENDIX C: FAST MEAN FIRST PASSAGE 
TIME CALCULATION 

In this subsection, we develop an approximation for 
computing the mean switching time very quickly. This 
approximation models the formation of a phase-slip train. 




bias current (fx A) 



FIG. 15: To compare approximate (thick solid lines) and exact 
(colored dots) methods for solving for the mean switching 
time delay-differential equation, we plot the mean switching 
rate T~^{Th, I) as a function of bias current I for various bath 
temperatures Tb = {L9, 1.7, 1.5, 1.3, 1.1, 0.9, 0.7, 0.5, 0.3} 
(from left to right) [l7|. The thin lines correspond to the 
phase-slip rate, and are shown for comparison. The parameter 
k in Eq. |C4j was set to 0.5. 



and is useful for fitting experimental data, where it is 
important to compute the mean switching time for a lot 
of points (Tb, /) quickly. 

In constructing this approximation we make several 
assumptions. First, we assume that the phase-slip trains 
are dilute in time, i.e., the wire spends most of its time 
at the temperature Tb, but very rarely there are trains of 
phase slips that heat up the wire. These trains are not 
overlapping, i.e., each train either leads to thermal run- 
away (a successful train) or the wire cools back down 
to Tb (an unsuccessful train). The train is considered 
successful if the wire temperature exceeds T* , as defined 
in Section ITVl 

In order to compute the switching rate, we must com- 
pute the probability for the formation of a successful 
phase-slip train S'(Tb, /) and multiply it by the phase-slip 
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rate r(Tb, /), which corresponds to the rate of formation 
of the first phase shp in a train. Thus the switching rate 
is given by 



T-i(Tb,/)-5(Tb,/)r(Tb,/). 



(CI) 



At this point, we make the additional assumption that 
the probabihty to form a successful train can be com- 
puted phenomenologically, as follows. Consider a phase 
slip in a wire that is at temperature T . Immediately af- 
ter a phase slip, the wire has temperature T -\- ri(T,I), 
but it is also cooling at the rate = (T — Tb) a(r, Tb). 
If the phase slip-train is to continue, there must be an- 
other phase slip within a time ^ r~^; otherwise the wire 
would cool to the bath temperature Tb and the phase-slip 
train would be unsuccessful. Applying this procedure to 
a chain of phase slips, we find that S'(Tb,/), the proba- 
bility to construct a successful phase-slip train, is given 
by 



N 



5(Tb,/)~n 



L (T, -Tb)a(T„Tb)" 



(C2) 



The mean switching rate, computed using the phe- 
nomenological model for the probability to form a suc- 
cessful phase slip train given by Eq. (|C2|1 . turns out to 
be too crude to give results that are quantitatively accu- 
rate, although, qualitatively, the exact results obtained 
by solving Eq. p2p are well reproduced. In order to im- 
prove accuracy, we take into account the fact that the 
phase-slip rate drops as the wire cools, and also introduce 
the tunable parameter fc, which characterizes how much 
the wire is allowed to cool before a phase-slip train is con- 
sidered to be unsuccessful. Consider a wire at tempera- 
ture Ti. In the absence of phase slips, we approximate 
the equation for the evolution of the wire temperature by 



dtT{t) = -{T-n)d,, 



(C3) 



where T(0) = T,, and = Ci(Ti,Tb). To parametrize 
the failure of a phase-slip train, we assume that a phase- 
slip train is unsuccessful if the temperature of the wire 
reaches the value 



T, 



fail 



(l-fc)T, + fcT,_i, 



(C4) 



where fc is a tunable parameter of order unity that should 
be chosen to minimize the difference between the exact 
[i.e., obtained from solutions of Eq. ((32|) ] and the phe- 
nomenological switching rates. Having defined T^^faii, we 
can define the time to reach it via 



ii.faii = -J- log 



T, ~ Tb 



(C5) 



Finally, taking into account the change in the phase-slip 
rate as the wire cools, as well as ii^taii, we modify Eq. (|C2[) 
to read 



N 



(C6) 

In Fig. [15] we compare the approximate switching rates 
obtained from Eq. (IClj) by using the phenomenological 
approximation Eq. (IC6p to the exact switching rate ob- 
tained from Eq. ((32| . We see that the phenomenological 
approximation is quantitatively very close to the exact 
switching rate. Therefore, to performing fits on exper- 
imental data we, in fact, use this phenomenological ap- 
proximation for the switching rate, as it can be computed 
much faster. 
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